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Many rotating machines such as compressors, turbines and pumps have long thin 
shafts with resulting vibration problems. They would benefit from additional damping 
near the center of the shaft. Magnetic dampers have the potential to be employed in 
these machines because they can operate in the working fluid environment unlike 
conventional bearings. This paper describes an experimental test rig which was set 
up with a long thin shaft and several masses to represent a flexible shaft machine. 

An active magnetic damper was placed in three locations: near the midspan, near one 
end disk, and close to the bearing. With typical control parameter settings, the 
midspan location reduced the first mode vibration 827., the disk location reduced it 
757. and the bearing location attained a 747. reduction. Magnetic damper stiffness and 
damping values used to obtain these reductions were only a few percent of the bearing 
stiffness and damping values. A theoretical model of both the rotor and the damper 
was developed and compared to the measured results. The agreement was good. 
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Equivalent stiffness 
Proportional gain 
Current stiffness 
Rate gain 

Controller total gain 
Position stiffness 
Mass 

Number of turns (per leg) 
Resistance 
Laplace variable 
Weight 

Vertical position 
Angular velocity 
Magnetic flux 
Permeability of free space 

Time constant 


INTRODUCTION 

Compressors, turbines and other rotating machines often have long thin shafts 
which give rise to high vibration problems. Currently, many compressors operate 
above the first [1] or second critical speed. Steam turbines also operate above 
several critical speeds. Seals inside compressors give rise to cross coupled 
stiffness terms which may drive a machine into large subsynchronous vibrations [2] . 
Nicholas et al. [3] report on the effects of flexible supports which tend to 
increase vibration levels in steam turbines. Generally it is desirable to increase 
the damping effects near the center of these long, thin shaft machines. One method 
of doing this is to place a damper inside the machine casing. 

The concept of using a damper inside these machines is not new. Such a damper 

would be located somewhere between the bearings and apply damping, and perhaps 
stiffness, to the rotor. What has been lacking until recent years is a damper which 
could operate in the environment of the working fluid. Conventional rolling element 
and fluid film bearings require an oil supply. A central damper in a compressor 
would have to be sealed so that the oil does not contaminate the gas. In most cases, 
the additional complications do not warrant the damper. In steam turbines, the 
temperature is so high that an oil damper could not be used inside the turbine. 

Magnetic dampers are now a real possibility for these types of machines. 
Commercial magnetic bearings [4,5] are increasingly being used for compressors and 
other machines. These do not require oil and can operate at high temperatures, thus 
removing the limitations indicated above. To date, the authors are not aware of the 
use of a magnetic damper in a field application in North America. 

A magnetic bearing can be used as a magnetic damper simply by having it carry no 

load. Thus the machine is completely supported in regular bearings. Typically these 
are oil lubricated bearings. If the damper is not active, then the machine reverts 
to its normal behavior without a damper effect. Some early work on simple rotors 
with active control was reported in [6,7] by Allaire et al. 
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Several studies have been carried out by Holmes et al. [8,9] on the suppression 
of vibrations in long transmission shafts. A magnetic damper was placed at various 
locations along the shaft to determine the reduction. These studies differ from the 
present paper because there were no large masses located along the shaft such as 
occur in compressors and turbines. 

A magnetic bearing/damper has been developed at the University of Virginia. It 
was used as the damper in this work, so the relevant publications are described in 
some detail here. The basic bearing geometry and control algorithm was described by 
Humphris et al. [10]. The effects of the bearing employed to support a flexible 
rotor were measured and reported in Allaire et al. [11-13]. Adjustment of the 
bearing stiffness properties via the control system gains allowed movement of the 
rotor first critical speed from approximately 800 rpm to 2,000 rpm, which is more 
than a factor of two. Variation of the bearing damping properties gave a reduction 
in rotor vibrations of as much as an order of magnitude. The first results using the 
magnetic bearing as a damper (not supporting any steady state load) were reported in 
[14] . Also, initial results of employing digital control to operate magnetic 
bearings were published in [15] . 

The purpose of this paper is to report on the use of an active magnetic damper 
in a multimass flexible rotor. Three large disks were placed on a shaft, supported 
in conventional bearings, to represent a multimass compressor or steam turbine. The 
magnetic damper was placed in various locations to observe the vibration reduction 
capabilities. The magnetic damper was placed near the midspan, near one end disk, 
and as close to the outboard bearing as possible. 


ACTIVE MAGNETIC DAMPER 

An active magnetic damper is essentially an active magnetic bearing with no load 
support capacity. The magnetic damper used in this analysis consists of four 
electromagnets located radially around the rotor shaft. The electromagnet currents, 
which determine the force in each magnet, were determined by a controller which 
maintained the shaft in the center of the damper by continuously monitoring the shaft 
position and adjusting the electromagnet currents accordingly. There were air gaps 
between the shaft and the magnets and the shaft did not contact the magnets. The 
dynamic properties of the damper (i.e. damping and stiffness) were electronically 
regulated by the controller. 

A diagram of the four electromagnets which made up the damper is shown in Figure 
1. They were located radially around a laminated ferromagnetic disk which was 
attached to the rotor shaft. The purpose of this laminated disk was to provide a 
continuous flux path between the magnet pole faces with a minimum of eddy current 
losses. 

Each electromagnet in the damper consisted of a solid soft iron core forming a 
horseshoe, with two pole faces cut to a diameter of 60.5 mm (2.38 inches). This gave 
a nominal radial clearance of 1.0 mm (0.040 inches). Each leg of the coil was wound 
with 920 turns of wire. All electromagnets were tne same with equal steady state 
currents. Presented in appendix A is the theory of the damper characteristics. 

Two proximitor (eddy current) probes located vertically and horizontally near 
the damper (unless otherwise specified) were used as rotor position sensors for input 
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to the electronic controller. Various tests were previously conducted to insure that 
the magnetic fields from the damper did not affect the readings of the probes. 


TEST RIG 

A four mass laboratory rotor was constructed with a 19.0 mm (0.75 inch) shaft 
supported in conventional sleeve bearings with a bearing span of 660.4 mm (26.0 
inches) . Three steel disks were placed four inches apart at the center of the shaft 
leaving 190.5 mm (7.5 inches) between each outboard disk and the closest support 
bearing. Each of these three disks was 139.7 mm (5.5 inches) in diameter, 25.4 mm 
(1.0 inch) in thickness, and weighed 28.48 N (6.4 lbs). The fourth mass was the 
laminated ferromagnetic damper disk which weighed 4.90 N (1.1 lbs) and was 58.4 mm 
(2.3 in) in diameter with a thickness of 25.4 mm (1.0 inch). The total rotor length 
and weight, including the laminated magnetic disk, was 762 mm (30.0 inches) and 107.2 
N (24.1 lb), respectively. 

This damper disk was placed at three different locations on the rotor, 
illustrated in Figure 2, depending on the desired location of the magnetic damper. 

The three locations were at the midspan near the center disk, near the third disk at 
the quarterspan of the rotor, and at a location just inboard of the outer bearing. 
These locations were chosen to demonstrate the effect of a damper at optimum and less 
than optimum rotor locations for control of one or more modes of vibration. These 
locations represent situations where there may be design constraints on an actual 
machine as to the placement of a magnetic damper on the rotor. As shown in Figure 2, 
shaft displacement sensors were placed at 3 different locations to monitor shaft 
motion. These locations will be referred to from left to right as the disk 1 probe, 
midspan probe and disk 3 probe, respectively. The four mass rotor will now be 
referred to as the 3- mass rotor, due to the relatively small mass of the laminated 
ferromagnetic disk compared to the three large steel disk masses. 

The 3- mass rotor was attached to an aluminum baseplate which was fixed to a 
large concrete block by anchor bolts. Shim stock was placed underneath the aluminum 
baseplate to correct a slight bow in the baseplate. The shaft was also twisted 
(bowed) with a significant shaft bow which could not be corrected. For example, when 
the shaft bow at the midspan probe was reduced to below .0254 mm (1 mil), the bow 
near the first disk was still close to .0508 mm (2 mils). 

The rotor was driven by a 1/2 hp electric motor which was rated at 10,000 RPM. 
The motor transmitted power through a flexible rubber coupling as shown in Figure 2. 

A total of three vertical and three horizontal proximitor probes were used to 
monitor shaft displacements. A separate proximitor probe was used as a key phasor. 
Signals from the probes were run to two digital vector filters, where amplitude and 
phase information at the probes was determined. One set of proximitor probes was 
also used as input sensors for the magnetic damper controller. These sensors were 
always located close to the damper. 

The three- mass experimental rotor was used to study the effects of a magnetic 
damper located at the three locations on the shaft. The magnetic damper, used in 
conjunction with a controller, exhibited some nonlinear effects above 4,000 RPM. 
Therefore this paper will be limited to operation below 4000 RPM (66.7 Hz). The 
first rotor bending mode is located at about 1600 RPM (26.7 Hz), with all other modes 
located at frequencies greater than 4000 RPM (66.7 Hz). 
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DAMPER AT MIDSPAN 


The first location of the magnetic damper to be investigated was the midspan 
location as shown in Figure 2a. This was expected to be the optimum position for 
control of first mode vibrations. 

The effect of increasing rate gain (K r ) or damping and proportional gain (Kg) or 

stiffness of the electromagnetic damper was studied at each location. The rate gain 
(K r ) was varied from a nominal value of 2.0 to a maximum value of 20.0. Due to 

stability problems, it was not possible to use this damper with pure damping or rate 
gain (K r ),and therefore, the controller was set with a nominal proportional gain (Kg) 

of 2.5, while the rate gain (K r ) was varied for analysis. Similarly, the 

electromagnetic damper could not be used with pure stiffness or proportional gain 
(K ), and the controller was set with a nominal rate gain (K ) of 5.0 while the 

O * 

proportional gain (Kg) was varied from 2.5 to 10.0. Finally, a comparison of the 
effect of the damper location for fixed values of rate gain (K r ) and proportional 
gain (K ) was carried out. 

o 

The effect of using the damper at the midspan is shown in Figure 3. The rotor 
was run with the damper off and with rate gain (K r ) equal to 2.0, 5.0 and 10.0. 

Proportional gain was held constant at 2.5. There was an 827. reduction in amplitude 
with a rate gain of 10.0 over the case with the damper off. 

An unbalance response model of the rotor was developed to compare to the 
measured data. Stiffness and damping coefficients were evaluated for the damper as 
described in Appendix B. Also, the unbalance level in the rotor was determined by 
balancing the rotor to a low level and then adding a known unbalance weight in the 
center disk. 

The forced response predictions for the magnetic damper at the midspan of the 
3- mass rotor are plotted in Figure 4 for cases with the magnetic damper off, and with 
the damper on with a rate gain (K r ) equal to 2.0, 5.0, and 10.0, respectively. A 

maximum reduction of 887. in the amplitude of vibration occurs theoretically with 
K r = 10 as compared to the value of 827. for the measured case. 

It may be seen from these results that even a small amount of damping can 
greatly reduce the vibration level. For K f = 2.0, the estimated damping value was 31 

N- sec/m (0.18 lb- sec/in), which produces a 497. decrease in amplitude at the rotor 
center. The estimated Searing damping was 1050 N- sec/m (6.0 lb- sec/in). Thus the 
actual midspan damping is only approximately 37. of the bearing damping. At K f = 10, 

the estimated damping was 122 N- sec/m (0.70 lb- sec/in) or approximately 127. of the 
bearing damping. This produced an 827. reduction in amplitudes of vibration. 

Figure 5 shows the experimental values for the midspan damper location as the 
proportional gain (K^) had the values 2.5, 5.0 and 7.5. The rate gain was held 

constant at K r = 5.0. The critical speed increased from 1600 rpm (26.7 Hz) to 2570 
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rpm (42.8 Hz) with = 7.5. The value of = 7.5 corresponds to a damper stiffness 

of 210,000. N/m (1200 lb/in). This stiffness value is only about 107, of the bearing 
vertical stiffness. 

The calculated forced response is shown in Figure 6. The predicted increase in 
critical speed is 387, for = 7.5. This is somewhat lower than the 607, measured 

increase in critical speed. Also, the theory predicted a steady increase in 
amplitude with increasing K , whereas the measured values decreased and then 

O 

increased. Overall, however, the agreement is good. 


DAMPER AT DISK THREE 

The second location of the magnetic damper is the disk three location as shown 
in Figure 2. This was expected to be the next best position for control of the first 
mode vibrations. 

The effect of using the damper at disk three is shown in Figure 7. The rotor 
was run with the damper off and the rate gain (K r ) equal to 5.0, 10.0, and 20.0. 

Proportional gain was held constant at 2.5. There is a 757. reduction in amplitude 
with a rate gain of 20.0 over the case with the damper off. 

An unbalance response model of the rotor with the damper at disk three was also 
developed for comparison with the measured data. The forced response predictions for 
this configuration are plotted in Figure 8 for cases with the magnetic damper off, 
and with the damper on, with a rate gain (K r ) equal to 5.0, 10.0, and 20.0, 

respectively. A maximum reduction of 857, in the amplitude of vibration occurred 
theoretically with K r = 20, as compared to the value of 757, for the measured case. 

Again, it may be seen from these results that even a small amount of damping can 
greatly reduce the vibration level. At K r = 20, the estimated damping was 175 

N- sec/m (1.0 lb- sec/in) or 177, of the bearing damping. This produced a 757. reduction 
in amplitudes of vibration. 

Experimental data was also taken for the disk three damper location with 
proportional gain (K ) values of 2.5, 7.5, and 10.0, while the rate gain (K r ) remains 

at 5.0. The critical speed increased from 1600 RPM (26.7 Hz) to 1920 RPM (32 Hz) 
with K =10.0. The value of K = 10.0 corresponded to a damper stiffness of 262,000 

(1500 fb/in). This stiffness value was only about 137. of the bearing vertical 
stiffness. 

The corresponding forced response calculations predicted a 257. increase in the 
critical speed frequency for K =10.0. This was somewhat higher than the measured 

O 

value of 197,. Once again the theory predicted a steady increase in amplitude with 

increasing K whereas the measured values decreased and then increased. In general, 

§ 

however, the agreement was good. 
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DAMPER AT BEARING 


The third location of the magnetic damper to be investigated was the bearing 
location as shown in Figure 2c. This was expected to be the least effective location 
for control of the first mode of vibrations. 


The effect of using the damper at the bearing was similar to the effect of the 
damper at other locations. However, experimental data was taken for the rotor when 
the rate gain (K r ) was equal to 5.0, 10.0, and 20.0, while the proportional gain (K^) 

equals 2.5. There was a significant 747, reduction in amplitude with a rate gain of 
20.0, over the case with the damper off whereas forced response calculations 
predicted only a 527, decrease. 


Similar data was taken when the proportional gain (K^) was set at 2.5, 7.5, and 

10.0, with the rate gain (K r ) equal to 5.0. There was a measured 117, increase in the 

first critical speed from 1620 RPM (27 Hz) with the damper off to 1800 RPM (30 Hz) 
with the proportional gain (K ) equal to 10.0. Forced response calculations 

o 

predicted a 67, increase in critical speed. 


EFFECT OF MAGNETIC DAMPER LOCATION 

A summary of the measured effect of magnetic damper location is shown in Figure 
9, when the proportional gain (K ) equals 2.5 and the rate gain (K ) equals 5.0. As 

noted, the optimum control of the first mode vibrations occurs when the damper is at 
the midspan location, with a 707, reduction. The next best location is at disk three, 
with a 617, reduction in vibration. While the least optimum location for control of 
the first mode of vibrations is the bearing location, there is still a 487, reduction 
in vibrations over the case with no damper. 

The corresponding forced response predictions are shown in Figure 10. There is 
a predicted 847, reduction when the damper is at the midspan and an 81% reduction when 
the damper is at disk three. Also there is a 347, predicted reduction when the damper 
is at the bearing. 

A summary of the overall measured responses at the first mode is shown 
graphically in Figure 11 with K = 2.5. This figure shows the effect of the damper 

o 

at various locations on the first mode response. It can also be noted that the 
vibration reduction approaches a minimum level as the damping is increased. 


CONCLUSIONS 

Overall, the magnetic damper achieves the desired result of a large reduction of 
first mode vibration at the center of the rotor. The largest reduction is 827, with 
the damper in the optimum location near the rotor center. These results indicate 
that even less than optimum placement or less than optimum control parameter setting 
values can produce very significant vibration reductions. This is very encouraging 
for the potential use of magnetic dampers in industrial applications. 
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It can also be noted, however, that there is a limit to the vibration reduction 
which can be attained by a single damper in the machine. Beyond a certain level, 
even a large increase in the stiffness and damping of the damper essentially produces 
very little further reduction of vibrations. 

Generally, the agreement between the theoretical modeling of the shaft/damper 
system compared to the measured results is good. The magnitudes are close and the 
trends in the results are consistent. This indicates that the design of magnetic 
dampers, done on a theoretical basis before construction, has a high probability of 
success. 


APPENDIX A 

ACTIVE MAGNETIC DAMPER 


This Appendix briefly discusses basic electromagnetic theory relating to the 
force and control of a magnetic damper. For an active magnetic damper there are four 
electromagnets located radially around the shaft. For the vertical direction there 
would be an identical top and bottom magnet and the force is 

.., 2 . 2 av 2 • 2 

H AN i, u AN i K 

-V- (1) 


In a damper, for steady state conditions, 


2 t ” 2 b 


The total force or load capacity in a magnetic damper is thus equal to zero. 


F total ■ » - 0 

This, of course, would not be true in a magnetic bearing. 

In an electromagnetic damper, two independent parameters can change -- position, 
x, and current, i. Define the position and current stiffness as 


AF 

Ky = - = Position Stiffness 


AF 

Kj = - = Current Stiffness 


The gap thickness and current expressions are 


h = h g - Ay 
i = i s + 

for small changes about the steady state. The total force expression, Equation (1), 
becomes 


F 


total 


/^(vAi) 2 


'// 0 AN 2 (i b -4i) 2 ~ 

. O s - ^) 2 . 


. (h g + Ay) 2 . 
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or approximating by the binomial expansion, 


total 


A 0 AN 2 i 2 
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For steady state conditions, 


i. = i. = i_ 

tbs 


so the total force equation reduces to 


total 


W~ 


4Av . 4Ai 
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The two stiffnesses K and K. can now be obtained. 


First setting Ai to zero, 


F = W - K y Ay = 


«AN 2 i 2 
r o s 


(r 1 ) 

S 


Recalling that W = zero for a damper, the total positive stiffness is 

(2) 


K = - 

y 

Second, setting Ay to zero, 


4/i AN 2 i 2 
r o s 


F = V - K.Ai = 


/* AN 2 i 2 
r o s 


(fi) 

S 


The total current stiffness is, 


K i = 


4/UN 2 i 
^0 s 

-n — 


( 3 ) 


Both the position and current stiffnesses for the combination top and bottom 
electromagnets are simply two times the position and current stiffnesses for a single 
electromagnet, respectively. 

The purpose of the electronic controller is to condition the position voltage 
signals from the proximitor probes and feed back the correct amount and phase of 
current to the magnets for optimum control of the rotor. The stiffness of the 
electromagnet damper is determined by the controller and can be varied by changing 
the proportional gain (K ) of the controller. The damping of the electromagnetic 

damper is also determined by the controller. The position voltage signal from the 
proximitor probes is electronically differentiated so that the damping can be 
determined by the amount of velocity or rate signal present. The damping can be 
varied by changing the rate gain (K ) of the controller. 
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A block diagram for the controller is shown in Figure 12. A feedback control 
diagram of the entire magnetic damper is shown in Figure 13. The transfer function 
for the electronics from the sensor amplifier through the position and current 
amplifier is 

G ( s ) = (U I - 1 s)(s+l'/T 3 )'(Ta , C5] 

Finally, the total transfer function for the magnetic damper is 

1 


(4) 


r ext 


Ms 2 + K 


1 + 


K/Cfif 

Ms 2 + K 


Simplifying 


'ext Ms 2 + K + K. G(s) 


(5) 


It is easily seen that adjustment of controller proportional and rate gains 
(K and K ) are used to adjust the stiffness and damping of the magnetic damper. 

An effective damper stiffness and damping can also be obtained. The transfer 
function for an equivalent mass- spring- damper system is 

JL_ = i 

ext Ms 2 * c eq s + K eq 

Equating this expression with Equation (5) gives 

K eq = K y + K i Re { G (j w )} 

-K i Im{G(jw)} (6) 

c eq “ w 

where s has been replaced by j w for synchronous vibrations. 


APPENDIX B 

DETERMINATION OF DAMPER DYNAMIC PROPERTIES 

A single mass rotor as shown in Figure 14 was used to experimentally evaluate 
the damping and stiffness characteristics of the magnetic damper for various settings 
of rate gain (K r ) and proportional gain (K ). The effective stiffness of the damper 

is determined by the proportional gain (K ) setting only. The effective damping of 

the magnetic damper is determined primarily by the rate gain (K r ) setting, but is 

also affected by the proportional gain (K ) setting. 

O 
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A rotor model of the single mass rotor was established, including support 
bearing characteristics, shaft bow, and unbalance in the rotor. The support bearings 
had a vertical stiffness of 4,030,000 N/m (23,000 lb/in) and a vertical damping of 
10,000 N- sec/m (57 lb-sec/in). This model was used in a forced response program and 
the predicted response was compared to experimentally measured rotor response. 

Various values of stiffness and damping were added to the rotor model at the 
electromagnetic damper location until predicted response (i.e. frequency and 
amplitude of the first mode response) matched the experimental rotor response to 
within approximately 57. . An example of experimental data and its corresponding or 
"matching" calculated response is shown in Figure 15 for the case when the 
proportional gain (K ) equals 2.5 and the rate gain (K r ) equals 2.0. In this manner 

stiffness and damping characteristics for the electromagnetic damper were determined 
for various values of rate gain (K r ) and proportional gain (K ) . The stiffness 

values are listed in Table 1 and the damping values are listed in Table 2. 

The damping and stiffness characteristics of the magnetic damper determined in 
this manner were used in the rotor model of the 3- mass rotor for forced response 
calculations. 
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Table 1 - Horizontal and Vertical Stiffness Values of the 
Magnetic Damper at Approximately 2500 RPM (41.7 Hz). 


Proportional Gain (K ) 

O 

0 

2.5 

5.0 

7.5 

10.0 

K x (N/n.) 

-35,000 

44,000 

190,000 

210,000 

260,000 

(lb/in) 

(-200) 

(250) 

(1100) 

(1200) 

(1500) 

K (N/m) 

-35,000 

70,000 

175,100 

210,100 

263,000 

(lb/in) 

(- 200) 

(400) 

(1000) 

(1200) 

(1500) 
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Table 2 - Horizontal and Vertical Damping Values (lb- sec/in) 
of the Magnetic Damper at 2500 RPM (41.7 Hz) 


N-sec, lb-sec. 


N-sec, lb-sec. 


N-sec, lb-sec. 


N-sec, lb-sec. 




Figure 1 Magnetic Damper Model 
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Figure 2 Schematic of 3-Mass Rotor with Magnetic Damper at Three Different Locations 
(Not to Scale) 
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AMPLITUDE - microns. 



Figure 3 Experimental Results for the Damper at Midspan (Increasing Rate Gain (K r )) 
at the Vertical Midspan Probe - = 2.5 



Figure 4 Forced Response Program Predictions for Damper at Midspan (Increasing Rate 
Gain (K )) at Vertical Midspan Probe - K = 2.5 
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SPEED, RPM x 10-3 


Figure 5 Experimental Results for the Damper at Midspan (Increasing Proportional 
Gain (Kg)) at the Vertical Midspan Probe - K r = 5.0 
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Figure 6 Forced Response Program Predictions for Damper at Midspan (Increasing 
Proportional Gain (Kg)) at Vertical Midspan Probe - K r = 5.0 
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Figure 7 Experimental Results for the Damper at Disk Three (Increasing Rate Gain 
(K )) at the Vertical Midspan Probe - K = 2.5 
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Figure 8 Forced Response Program Predictions for Damper at Disk Three (Increasing 
Rate Gain (K )) at Vertical Midspan Probe - K = 2.5 
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Figure 9 Experimental Results at Vertical Midspan for the Damper at Three Locations 
(Proportional Gain (K ) equal to 2.5, Rate Gain (K r ) equal to 5.0) 

O 
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Figure 11 Response Amplitude at First Mode vs. Rate Gain (Damper at Midspan, Disk 
Three, and Bearing Location) 
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Figure 12 Block Diagram of Control Circuit 
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MASS - MAGNET 



Figure 13 Diagram of Magnetic Damper System for Controller 



Figure 14 Schematic of Single Mass Rotor with Magnetic Damper (Not to Scale) 
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Figure 15 Comparison of Calculated and Experimentally Measured Forced Response 
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